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Quantum networks promise to provide the infrastructure for many disruptive 
applications, suchas efficient long-distance quantum communication and 
distributed quantum computing”. Central to these networks is the ability to 
distribute entanglement between distant nodes using photonic channels. Initially 
developed for quantum teleportation** and loophole-free tests of Bell’s inequality~*, 
recently, entanglement distribution has also been achieved over telecom fibres and 
analysed retrospectively”*. Yet, to fully use entanglement over long-distance 
quantum network links it is mandatory to knowit is available at the nodes before the 
entangled state decays. Here we demonstrate heralded entanglement between two 
independently trapped single rubidium atoms generated over fibre links witha 
length up to 33 km. For this, we generate atom—photon entanglement in two nodes 
located in buildings 400 m line-of-sight apart and to overcome high-attenuation 
losses in the fibres convert the photons to telecom wavelength using 
polarization-preserving quantum frequency conversion’. The long fibres guide the 
photons to a Bell-state measurement setup in which a successful photonic projection 
measurement heralds the entanglement of the atoms”. Our results show the 
feasibility of entanglement distribution over telecom fibre links useful, for example, 


for device-independent quantum key distribution 


13 and quantum repeater 


protocols. The presented work represents an important step towards the realization 
of large-scale quantum network links. 


Sharing entanglement between distant quantum systems is a crucial 
ingredient for the realization of future quantum networks. Photons 
are the tool of choice to mediate entanglement distribution, typically 
either by controlled light-matter interaction with local memories», 
or, as it also will be used here, by entanglement swapping from two 
pairs of entangled photon-memory states””’*"®, Innovative applica- 
tions of such networks include distributed quantum computing” and 
device-independent quantum key distribution” ©. As attenuation losses 
in the distribution process are inevitable, quantum repeaters will be 
essential to efficiently distribute entanglement by intermediate nodes. 

To minimize absorption along the quantum channel and thus to 
maximize the distance between neighbouring nodes in quantum 
networks using the readily available fibre infrastructure, it is neces- 
sary to operate at telecom wavelengths. Using quantum frequency 
conversion (QFC)*°™, light-matter entanglement distribution at the 
low loss telecom band has recently been demonstrated for various 
types of quantum memory, including NV-centres, ions, atoms and 
atomic ensembles”*’, even over tens of kilometres of fibre’”’. This 
was mainly enabled by new quantum frequency converters”>”*, which, 
while preserving the photonic polarization, have reached external 
device conversion efficiencies as high as 57% (ref. ’). 


For future quantum communication and repeater scenarios, it is vital 
that the nodes are independent and distant, use long-lived quantum 
memories and, at the same time, provide the availability of heralded 
entanglement, that is, there is a well defined signal available that the 
entanglement distribution succeeded. So far, this has been limited 
to fibre lengths up to 1.7 km (refs. °°). Recently, great progress was 
made by demonstrating entanglement between atomic ensembles’ and 
multimode solid-state quantum memories using telecom wavelength 
photons, however, this had limited memory storage times and did not 
use independent nodes. 

Here we report on the distribution of entanglement between two 
remote quantum nodes—*’Rb atoms trapped and manipulated inde- 
pendently at locations 400 m apart—generated over fibre links witha 
length of up to 33 km. The scheme begins with entangling the spin state 
of anatom with the polarization state of a photon in each node. Subse- 
quently, the photons emitted by the atoms at 780 nmare converted to 
telecom wavelengths and transferred over several kilometres of fibre to 
amiddle station, where a Bell-state measurement is performed to swap 
the entanglementtothe atoms. We analyse the heralded entanglement 
between the atoms for different fibre link lengths using correlation 
measurements of the atom-atom state along three bases. The atoms are 
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Fig. 1| Schematic of the experimental setup. In each node, located in 
buildings 400 mapart, asingle *’Rb atomis loaded in an optical dipole trap. 
Both atoms are synchronously excited to the state 5”P;/.|/F’=0, m;-=0)to 
generate atom-photon entanglement in the subsequent spontaneous decay. 
The single photons emitted at a wavelength of 780 nmare collected using 
high-numerical aperture objectives and coupled into single-mode fibres 
leading to the QFC devices. There, they are converted to telecom wavelength 
(A=1,517 nm) by difference frequency generation ina periodically poled 


analysed after a delay corresponding to the two-way communication 
time to the middle station over the full fibre link length to realistically 
evaluate the performance for long fibre lengths. 


Quantum network link 


Our experiment consists of two similar, independent nodes anda mid- 
dle station, all of which are located in different laboratories, as illus- 
trated in Fig. 1. The shortest possible fibre connection from node1 
(node 2) to the middle station equals 50 m (750 m), longer fibre links 
are realized by inserting extra fibres on spools. The fibre length to the 
middle station is denoted as L, (L,), with the totallinklengthZ =L,+L,. 
In bothnodes, asingle, optically trapped *’Rb atom acts as aquantum 
memory”, in which a qubit is encoded in the Zeeman substates of the 
5°S12|F=1, m-=+1) ground state, with m;=+1 and m,= 1 further 
denoted as|*), and|V),, respectively. 

The experimental sequence starts by generating atom-photon entan- 
glement in each node”. For this purpose, the atoms are prepared in the 
initial state 3 SylF= 1, m; = O) by optical pumping and excited to the 
state 5P3/.|F’ = 0, m,,= 0). During the spontaneous decay back to the 
ground state, the atomic spin state becomes entangled with the polar- 
ization of the respective emitted photon at 780 nm due tothe conserva- 
tion of angular momentum. This results in the entangled atom- 
photonstate |W) 4p = 1/2 (|v) IL) + |*),1R)) =1//2 (|v) eH) +1) IV), 
where |L) and |[R) denote left- and right-circular photonic polarization 
states, |H) and |V) denote horizontal and vertical linear polarizations. 
Acustom made high-numerical aperture objective is used to collect 
the atomic fluorescence into a single-mode fibre with an efficiency of 
1.0% (1.1%) after an excitation attempt in node 1 (node 2). 

Photons witha wavelength of 780 nm now would suffer an attenu- 
ation bya factor of ten after propagation through 2.5 km of fibre. To 
overcome such high-attenuation loss, we use polarization-perserving 
QFC to transform the wavelength of the photons to the telecom S 
band, in which one expects attenuation by a factor of ten after only 
about 50 km transmission. The QFC is realized by mixing the pho- 
tons with a strong pump field at 1,607 nm inside a non-linear wave- 
guide crystal, converting the wavelength to 1,517 nm by difference 
frequency generation. Various spectral filtering stages, including 
a narrow band filter cavity (27 MHz full-width at half-maximum), 
separate the converted single photons from the strong pump field 
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lithium niobate (PPLN) waveguide located in a Sagnac interferometer type 
setup, sucha configuration fully maintains the polarization quantum state of 
the photon. The converted photons are guided toa middle station by fibre links 
with lengths up to 16.5 km, where the entanglement is swapped to the atoms by 
aBSM. After successfully generating atom-atom entanglement, the atoms are 
analysed independently by areadout pulse of which the polarization, set bya 
half-wave plate (HWP) and quarter-wave plate (QWP), defines the measurement 
setting. PC, polarization controllers. 


and the anti-Stokes Raman background originating from this field. 
In the shortest fibre configuration, this results in a background of 
roughly 160 and 170 cps registered at the middle station for light 
from nodes 1 and 2, respectively. Both converters achieve an external 
device efficiency of 57%. The pump light is conveniently distributed 
to the nodes using the telecom fibre network and hence ensures 
indistinguishable frequencies of the single photons after conver- 
sion. For more details about the QFC system and an analysis of the 
atom-photon entanglement distribution at telecom wavelengths, 
see ref. ° and Methods. 

After frequency conversion, the photons are guided to the middle sta- 
tion with fibres of different lengths where an interferometric Bell-state 
measurement (BSM) swaps the entanglement to the atoms°”’. The 
fidelity of the BSM, and hence of the entanglement swapping protocol, 
is determined by the photons temporal, spectral and spatial indistin- 
guishability*. This is optimized by different means. First, the photons 
impinge ona balanced, single-mode fibre beam splitter to guarantee a 
perfect spatial overlap. Second, the entanglement generation process 
in the nodes is synchronized to <300 ps, which is much smaller than 
the coherence time of the photons determined by the lifetime of the 
excited state of 26.2 ns. And third, polarization drifts in the long fibres 
are compensated using an automated polarization control*. The pho- 
tons are detected with four superconducting nanowire single-photon 
detectors (SNSPDs), which all have an efficiency of >85% and a dark- 
count rate of <65 cps. 

The BSM setup used here analyses the photons in the H/V basis and 
hence heralds the following two Bell states of the atoms 


V2) =U) )e# I) eI) (1) 


Two-photon coincidences are triggered within a hard-wired 208 ns 
long window that sends a heralding signal back to the nodes. The signal 
is delayed electronically by t= ¢/=cto simulate the signalling time back 
to the nodes, where Zc approximates the speed of light in an optical 
fibre and? equals L, or L, for node1or 2, respectively. Although lower- 
ing the observed final fidelity, this delay is essential to study the per- 
formance of the quantum network link ina realistic scenario. 

The quantum state of the atomic memories is finally analysed with 
a state-selective ionization scheme, in which the state selectivity is 
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Fig. 2 | Characterization of the atom-atom entanglement for afibre length 
ofL = 6 km.a, Detection time histogram of the photons originating from 
node land 2 relative to the time of excitation in node 1. For the indicated 
acceptance window, we observe a SBR of 58 and 65 for nodes 1and 2, 
respectively. The temporal overlap of the two photons is >0.98. b, Two-photon 
interference based on the Hong-Ou-Mandeleffect. Shown is the normalized 
(norm.) number of wrong coincidences for various time differences between 


controlled by the polarization of a readout laser pulse’. Both memo- 
ries have a coherence time T, of roughly 330 ps, which is achieved 
by active compensation of external magnetic fields (<O.5 mG) and 
applying a bias field of tens of milligauss along the y axis. Currently, 
the coherence time is limited by magnetic field fluctuations along the 
bias field direction and a position-dependent dephasing originating 
from longitudinal field components of the strongly focussed dipole 
trap. For more details and a simulation of the limiting decoherence 
mechanisms, see Methods. 


Entangling atoms using telecom photons 


The evaluation of the entanglement distribution over long fibre links 
is detailed first for a fibre configuration of L = 6 km. The entanglement 
generation rate is determined by two factors: the probability of gen- 
erating an entangled state between the atoms after a synchronized 
excitation attempt, predominantly reduced by the photon collection 
efficiency in the nodes, amounts to 7 = 3.66 x 10° and the repetition 
rate of R = 30.8 kHz, whichis mainly limited by the communication time 
between the nodes. This leads to an event rate r= 1/9 s “for coincidence 
detections in the full hard-wired window. With a duty cycle of the setup 
of roughly half, including the fraction of time that an atom is presentin 
both traps of 0.60 (5), we observed N = 10,290 entanglement generation 
events in 54 h. For more details about the entanglement generation 
rate and duty cycle, see Methods. 

Facing Raman background from two QFC devices in addition to 
detector dark counts, a 70 ns photon acceptance window is applied 
during the data evaluation, as shown in Fig. 2a. This results in a 
signal-to-background ratio (SBR) of 58 (65) for detecting a single pho- 
ton from node 1 (node 2), whichis substantially higher than in a previ- 
ously reported work (ref. ®) thanks to a more favourable pump-signal 
frequency combination with respect to the Raman background. For 
the coincidence detections this leads to a SBR of 48, while accepting 
roughly 65% of the recorded events. 
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the two photon wave-packets (Ar). c, Atom-atom state correlations for the|W" ) 
(top) and|*) (bottom) states. The correlation probability of the measurement 
outcome in the nodes is shown in blue, whereas the anticorrelation 
probabilities are marked in red. The dataare fitted with sinusoidal functions 
resulting in an estimated fidelity for|W )and|W*) relative toa maximally 
entangled state of 0.826(18) and 0.806(20), respectively. The error bars in this 
and all subsequent figures indicate statistical errors of one standard deviation. 


The quality of the two-photon interference of the converted pho- 
tons is quantified by the relative occurrence of wrong detector coin- 
cidences”. These coincidences, that is, (V,,V,) and (H,,H,), should not 
occur for perfectly interfering, fully unpolarized photons. For tempo- 
rally well overlapping photons (At = 0), but without background correc- 
tion, this results in an interference contrast of 0.955(7) (Fig. 2b). Here, 
two photon emission events from one atom reduce the indistinguish- 
ability of the interfering photons by changing the temporal photon 
shape™. By rejecting early detection events, this effect is reduced at 
the cost of alower event rate. For more details, see Methods. 

To evaluate the atom-atom entanglement we measured the atomic 
spin states in the two linear bases, X and Y. For this, the polarization 
analysis angle in node 2 was set to 6 = 0° = Xand B= 45° = Y, whereas 
the analysis angle in node 1 was varied over 90° in steps of 22.5° starting 
froma=0°=Xanda=45° = Y, respectively. The atom in node 1 (node 2) 
was analysed at ¢, = 28.5 ps (t, = 35.5 ps) after the respective exci- 
tation pulse. The resulting atomic state correlation probabilities 
Poorer = (N28 + N%-P)/N*8 and anticorrelation probabilities P,.o,.= 
(NP +N&P)/N“8 are shown in Fig. 2c. The data are fitted with 
sinusoidal curves giving average visibilities of V = 0.804(20) for |W) 
and V = 0.784(23) for |W"). To estimate the state fidelity we need to 
consider that the third ground state(575, ./F = 1, mp = 0)) can be popu- 
lated and hence the atomic states evolve effectively in a3 x 3 state space. 
Therefore, a lower bound on the fidelity is given by F= ; + SV 
(Methods), resulting in fidelities of F>0.826(18) for |W’) and 
0.806(20) for|*), relative to a maximally entangled state. (When also 
accepting the remaining 35% of the coincidences in the hard-wired 
acceptance window, the average observed fidelity for the states |W*) 
equals 0.772(10), with a SBR of 22, and an interference contrast of 
0.895(8).) The analysis angles chosen allow to evaluate the Clauser— 
Horne-Shimony-Holt S value® for the settings a = 22.5°, B= 0°; 
a’ = 67.5°, B= 0°; a’ = 67.5°, B= 45° and a” = 112.5°, 6’ = 45°, whereby a” 
replaces a = 22.5°. This results in an observed value of S = 2.244(63), 
violating the classical limit of two with 3.90. 
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Fig. 3 | Observation of heralded atom-atom entanglement generated over 
long fibre links. a,b, Correlations between the measurement results when 
analysing the generated atom-atom Bell states|W” ) (a) and|W*) (b) for various 
link lengths. For each link, the states were analysed in three conjugate bases 
(node 1, node 2), whereby the correlation probability of the measurement 
result in the nodes equals P.o,-= (Nr. + Nyy )/(Npn t+ Nay + Ny. + Nyy). For the 
different fibre lengths (short to long) 4,281, 4,271, 4,153 and 3,022 
entanglement events were recorded within measurement times of 11, 65, 97 
and 175 h. Of these events, 62-72% were within the two-photon coincidence 
acceptance window, resulting in N= 185 to 225 events per data point. 
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To determine the quality of entanglement for long fibre links, we per- 
formed a series of measurements generating and observing atom-atom 
entanglement in fibre configurations with a length L = 6, 11, 23 and 
33 km. For longer links, the event rate reduces because of both the 
signal attenuation around 0.2 dB km“ and the longer communica- 
tion times (Methods). At 33 km, this results in a success probability 
of 7 =1.22 x 10°, a repetition rate of R = 9.7 kHz and an event rate of 
r=1/85s". The entanglement fidelity relative to maximally entangled 
states was analysed by measurements along three bases (X, Y, and Z). 
Fig. 3 shows the probability of correlated measurement outcomes 
in the nodes for each measurement setting combination and fibre 
configuration. 

The fidelity of the observed states is estimated by first determining 
the contrast in the three measurement bases independently. This is 
done by taking the absolute difference of the two measured correlation 
probabilities, F, = |P, ,—P-, ,|for k € {X, Y, Z}, from which the average 
contrast is computed as EF = (Ey + Fy + Ez)/3. When averaging over the 
observed states |W*)—which show in our measurement precision sim- 
ilar visibilities—and again assuming the 3 x 3 state space, this results 
in alower bound on the fidelities F= 0.830(10), 0.799(11), 0.719(12) 
and 0.622(15) relative to maximally entangled states for L equals 6, 11, 
23 and 33 km, respectively. The estimated fidelity for the 6 km fibre 
configuration is in good agreement with the fidelity estimated from 
the measurements in two bases presented before (Fig. 2c). Moreover, 
allobserved fidelities clearly exceed the bound of 0.5 and hence clearly 
witness an entangled state. 

The observed fidelities are shown in the dependence of the fibre link 
length in Fig. 4, where the measured atom-photon entanglement 
fidelities for states shared between the nodes and the middle station 
are also shown for comparison (Methods). The atom-atom state 
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Fig. 4| Observed entanglement fidelity for various link lengths. Overview 
of the observed atom-atom fidelities for different fibre configurations (blue). 
For completeness, the observed atom-photon fidelities of the states between 
the nodes and the middle station (red), with = 2L,, are given (Methods). The 
solid lines are simulations based ona model taking into account the 
decoherence of the atomic memories (Methods). The grey points are reference 
measurements of the atom-atom state decoherence without long fibres but 
with corresponding readout delay. 


fidelity for different fibre lengths A(L) is modelled on the basis of 
simulations of the generation and evolution of the two atom-photon 
states, which are shown with solid lines. We estimate the visibility of 
the atom-atom state by the product of the two atom-photon visibili- 
ties and the interference contrast of the BSM*. Evidently, decoherence 
of the atomic states dominates the loss in fidelity for longer fibre links. 
For L = 33 km, the states were analysed 171 and 178 ps after excitation 
in nodes 1 and 2, respectively, which approaches the coherence time 
ofthe states. By contrast, the SBR is robust to an increase in fibre length 
because both the single photons as well as the QFC background are 
attenuated in the long links. A minor reduction in SBR (42 for L = 33 km) 
is explained by relatively more detector dark counts and can be solved 
by installing detectors with lower dark counts. Also, polarization drifts 
are comparably well compensated in all link configurations (Methods). 

To verify that the memory decoherence limits the loss in fidelity for 
long fibre links, we performed a series of measurements without further 
fibres inserted, however, withthe memory readout times electronically 
delayed according to the long fibre links. The observed fidelities are 
shown in grey in Fig. 4 at L = Zc(t,+ t,)/2 (matching the two-way com- 
munication time to the middle station for distance L) and show, within 
the measurement accuracy, no difference in observed fidelity com- 
pared to the configuration with long fibres. 


Discussion and outlook 


The results clearly indicate the feasibility of turning to large-scale quan- 
tumnetworks and of increasing the line-of-sight separation of the nodes 
to tens of kilometres by using efficient QFC. To evaluate the perfor- 
mance of entanglement generation in future quantum networks the so 
called quantum link efficiency was introduced recently”. It is defined 
as the ratio of the entanglement generation rate over the entanglement 
decay rate and describes howefficiently one can use entanglement asa 
resource in future quantum networks. Ideally, it should exceed one, that 
is, entanglement is available on demand as it is generated faster than 
it decays. However, the link efficiency decays rapidly with length due 
to both the exponential decrease of the signal detection probability, 


but, even more dominantly for the link lengths realized here, also due 
to the waiting times for classical communication between the nodes 
(Methods). In this proof-of-principle demonstration, the link efficiency 
was at the 10° to 10° level, mainly due to the low photon collection 
efficiency, relatively short coherence times and long fibre links. 

To raise the link efficiency, several improvements are feasible: 
the coherence time of the atomic memories can be increased by 
the implementation of, first, a new trap geometry to mitigate the 
position-dependent dephasing and second, a state-transfer to a qubit 
encoding less sensitive to magnetic fields** making possible coherence 
times in excess of 5 ms, and thus distances of 100 km without notable 
memory decoherence. On the entanglement generation side, opti- 
cal cavities could enhance the fluorescence collection efficiency into 
fibres. For single rubidium atoms, photon collection efficiencies up to 
11% have been demonstrated”. Yet, even for an almost ideal experimen- 
tal platform the entangling rate cannot be increased arbitrarily high. 
To further increase the event rate on the long run it will be mandatory 
to parallelize the entanglement distribution to regain the scalability. 
For the neutral atom platform this could be achieved by realizing atom 
trap arrays” and hence potentially increase the entanglement genera- 
tion rate by orders of magnitude. This concept will also lay the basis to 
realize quantum repeater network nodes performing entanglement 
purification and BSMs by local atom-atom gate operations inthe array, 
for instance, using Rydberg interactions”. 

To conclude, using efficient telecom interfaces in our nodes enabled 
the generation of heralded entanglement between two atomic quantum 
memories over fibre links with a length as long as 33 km. The results 
clearly show that improvements on the memory coherence time and 
entanglement generation rate are mandatory, but will allowto entangle 
two atomic quantum memories with a fidelity better than 80% over 
fibres links up to 100 km, thereby paving the way towards long-distance 
entanglement distribution for future quantum networks. 
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Methods 


Entanglement generation and single atom state readout 
schemes 

The atom-photon entanglement generation process is visualized in 
Extended Data Fig. la. It starts by preparing the atom in the state 
3 SyalF= 1,m,;=0) by optical pumping with an efficiency of 80%. 
Next, a short laser pulse with 21 ns full-width at half-maximum 
excites the atom to the state 5”P,/,|F’= 0, m;, = 0). In the subsequent 
spontaneous decay, the atomic spin state becomes entangled 
with the polarization state of the emitted photon, denoted as 
IW) 4p = 1/2 (IW) IL) + 1) IR) = 1/2 (IY) lH) + [*) x 1V)) , where |H)= 
i/.J2 (IL) ~ IR)),|V)=1/J2 (IL) + IR)), |v) y = i/V2 (IP) — Y)z), ANA |) = 
1/J2 (It), +|¥)z). 

The atomic state readout process is visualized in Extended Data 
Fig. 1b. After a successful atom-atom entanglement generation event, 
the atomic spin states of the atoms are individually analysed using a 
state-selective ionization scheme. This scheme starts by transferring 
a selected atomic qubit state superposition from the ground state 
5°S,/2|F = 1)to the excited state 5”P, .|F’ = 1) with light at 795 nm (‘read- 
out light’). Simultaneously, the excited state is ionized using a bright 
laser pulse at 473 nm. If the atom possibly decays to the state 3 SylF= 2) 
before ionization, as indicated with the right-most grey arrow in 
Extended Data Fig. 1b, a780-nm cycling pulse will transfer it to the state 
5?P3:|F’ = 3), which is also ionized. The state readout is completed by 
fluorescence collection on aclosed atomic transition to check whether 
the atom is still present in the trap or not. The fidelity of atomic state 
readout operation is 96%. 

The measurement basis of the atomic qubit state is controlled by 
the polarization of the readout light pulse, which is defined as 
x= cos(a)V+e  sin(a)H. Accordingly, two orthogonal atomic qubit 
state superpositions can be derived of which one is transferred to the 
excited state by the readout pulse (bright-state) and the other is not 
(dark state), given as 


I )pright-state = €0S(@) (IV), —1),) + sin(aer# (It), +%),) (2) 


J2 


I) bork-state = sin(a) (WY): -*),)- costae (WV); +t). @) 


Note that all states except the dark state are excited and hence 
ionized, for example, population in the state 5’5, /F=1, m, =0) is 
always ionized. This makes the readout scheme a projection measure- 
ment onto the dark state. 

Anintuitive example of the state selectivity is the case of ao*-polarized 
readout pulse. In the z basis, as shown in Extended Data Fig. 1b, this 
pulse will excite an atomin the state|v ),tothestate 5°P,/.|F’ = 1, mp, = 0), 
however, an atom in the state|*), will not be excited because of the 
absence of the state 5”P,.|F’=1, mp, =2). 


Atom-photon entanglement distribution at telecom 
wavelength 

The polarization-preserving QFC devices used in this work are described 
in detail in refs. °*". In contrast to previous work, herea more favourable 
pump-signal frequency combination is selected with respect to the 
Raman background: 1,607-1,517 nm instead of 1,600-1,522 nm. This 
increases the SBR by a factor of four and allows to install a QFC device 
in both nodes without being limited by the SBR. 

Because the quality of the entanglement shared between the two 
nodes directly depends onthe fidelity of the two entangled atom-pho- 
ton pairs, we individually characterize the atom-photon entanglement 
generated in both nodes. The generated states are analysed using the 
same fibre configurations and atomic readout times as during the atom- 
atom entanglement measurements presented in the main text. For an 


overview, see Extended Data Table 1. Note that ahigh fidelity atomic state 
readout can only be made after a full oscillation period of the atom inthe 
dipole trap, which equals 14.3 and17.8 ps for nodes 1and 2, respectively. 

The atom-photon entanglement fidelity is analysed following the 
methods in ref. °, in which the atomic readout time is now delayed to 
allow for two-way communication to the middle station for each node 
over the respective fibre length. The polarization of the photons are 
measured in two bases, H/V (horizontal/vertical) and D/A (diagonal/ 
antidiagonal), that is, X and Y, whereas the atomic analysis angle was 
rotated over angles including these bases. The atom-photon state 
correlations are shown for node 1 in Extended Data Fig. 2 and for node 
2 in Extended Data Fig. 3. 

For the fibre configuration L = 6 km, thatis, L, = 2.6 kmandL, = 3.3 km, 
we find atom-photon state fidelities of 0.941(5) for node land 0.911(6) 
for node 2, relative to amaximally entangled state, which are mainly lim- 
ited by the atomic state readout and entanglement generation fidelity. 
For longer fibre lengths, the fidelity of the entangled state reduces 
due to magnetic field fluctuations along the bias field direction and 
the position-dependent dephasing. 


Modelling of the quantum memory decoherence 

Inbothnodes, asingle ®’Rb atom is stored in an optical dipole trap where 
a qubitis encodedinto the states SS 2\F= 1,m,=+1). The dipoletrapis 
operated at A op; = 850nm with typical trap parameters of, for example, 
for node 1, atrap depth U, = 2.32 mK and beam waist w, = 2.05 pm. The 
qubit evolves effectively inaspin-1systemas the state 5°5, iF =1, mp = 0) 
could also be populated. The state fidelity is influenced by two factors: 
the first one is the a.c. Stark shift originating from the dipole trap and, 
second, the Zeeman effect arising from magnetic fields. 

To model the dephasing of the quantum memories, we simulate the 
evolution of this spin-1 system while the atom is oscillating in the dipole 
trap, affected by longitudinal polarization components and external 
magnetic fields*. For this, we first randomly draw a starting position 
and velocity of an atom from a 3D harmonic oscillator distribution in 
thermal equilibrium. Second, the motion of the atom is simulated ina 
realistic Gaussian potential resulting in an atomic trajectory for which 
the evolution of the atomic state is calculated on the basis of the local 
optically induced and external magnetic fields. Finally, this is repeated 
for alarge number of trajectories in which the averaged projection for 
all trajectories yields the simulation result. 

The model takes the following independently measured inputs: (1) 
the trap geometry specified by the beam waist @, which is obtained 
from knife-edge measurements of the dipole trap beam focus in two 
dimensions**; (2) the trap depth U,, determined by measurements 
of the transverse trap frequency using parametric heating** and the 
atomic state rephasing period” and (3) the atomic temperature 7, 
modelled as a Boltzmann distribution that is measured by the release 
and recapture technique®. Inputs 1 and 2 define the position, amplitude 
and phase of the longitudinal polarization components, whereas inputs 
1-3 characterize the atomic trajectories. Furthermore, we include a 
uniform magnetic field along three directions with shot-to-shot noise 
following Gaussian distributions. 

Extended Data Fig. 4 shows simulation results and measurement data 
of the state evolution in node 1 for varying state readout orientation 
and time. The model accurately predicts the evolution of the measured 
atomic states and shows that the memory storage time is limited by 
magnetic field fluctuations onthe order of <O.5 mG along the bias field 
direction in addition to the position-dependent dephasing due to the 
longitudinal field components of the strongly focussed dipole trap. The 
simulation results presented in the main text consider the envelope of 
the found oscillating state evolution in three bases. 


Experimental sequence 
The entanglement generation sequence is visualized in Extended 
Data Fig. 5. The sequence starts by trapping an atom in both nodes. 


For this, a single atom is loaded from a magneto-optical trap into a 
tightly focussed dipole trap, which takes roughly 1s (2 s) fornode1(2). 
Every entanglement generation try consist of 3 p's optical pumping 
(80% efficiency) and an excitation pulse (Gaussian laser pulse with a 
full-width at half-maximum of 21 ns) to generate atom-photon entan- 
glement in the following decay. Subsequent to each try, a waiting time 
is implemented to cover the propagation time of the photons in the 
long fibres. After 40 unsuccessful tries, the atoms are cooled for 350 ps 
using polarization gradient cooling. The lifetime of the atoms inthe trap 
during the entanglement generation process is 4 s (6s) for node 1 (2). 

To verify whether bothtraps still store a single atom during the entan- 
glement generation tries, the process is interrupted after 200 ms to 
check the presence of the atoms. For this, a microelectromechanical 
systems, a fibre-optic switch is installed in each node at the SM-fibre that 
is used for the atomic fluorescence collection from the atom trap. The 
switches guide the atomic fluorescence either to the QFC devices during 
the entanglement generation tries, or to an avalanche photodiode 
located at each node during 40 ms of fluorescence collection. Note 
that the SNSPDs of the BSM cannot be used for this purpose because 
they are behind narrowband spectral filters. 

Regular maintenance tasks lower the duty cycle of the experiment 
to roughly half for all link lengths. This includes the fraction of time 
required to simultaneously load an atom in the traps (0.40(5)), the 
fraction of time used to verify if both traps still store a single atom 
during the entanglement generation tries (0.18), and the fraction of 
time used to compensate polarization drifts of the long fibres (0.05). 


Entanglement generation rate 
The atom-atom entanglement generation rate ris given by 


r=nR, (4) 


here 7 equals the success probability for each entanglement generation 
try and Ris the repetition rate of the entanglement generation tries. 
Both 7 and R are dependent on the link length L. In the following, we 
assume a two-node setup witha middle station halfway the nodes, that 
is,L=L,+L,withL, =L,, where L, (L,) equals the link length from node1 
(node 2) to the middle station. 

The success probability of an entanglement generation try is given by 


n(L) =n, 10°" = O(exp(-L)), (5) 


here 7,-, denotes the success probability for a setup witha zerolength 
link and approximates 5.0 x 10° for the presented apparatus. This 
includes the photon collection efficiencies in both nodes after an 
excitation attempt (1.0 and 1.1%), the transmission of the microelec- 
tromechanical system switches (85%*), the efficiency of the frequency 
conversion devices (57%*), the single-photon transmission efficien- 
cies of the spectral filtering cavities (81%*), the fibre coupling to the 
single-photon detectors (90%*), the single-photon detector efficien- 
cies (85%*) and the fraction of distinguishable Bell states of 2/4 (*these 
efficiencies should be included twice). The attenuation rate in opti- 
cal fibres is denoted by ain units of dB km“, which is reduced using 
polarization-preserving QFC to telecom wavelength from 4.0 dB km™ 
at 780 nm to 0.2 dB km‘ at 1,517 nm. 
The repetition rate of the entanglement generation tries equals 


1 = 1 
T Tot $L/ (3c) 


R(L) =O(L"), (6) 


where Tis the period of an entanglement generation try. The period 
of an entanglement generation try for a link with L = 0 is denoted by 
T,-9 and equals 12 1s for the presented apparatus (R(O) = 8.3 x 10* s). 
This includes the initial state preparation (3 ps), the entanglement 
generation (200 ns), and the duration of the polarization gradient 


cooling per try (350 p1s/40) to counteract the introduced heating dur- 
ing the state preparation and entanglement generation tries. The sec- 
ond term inthe denominator gives the communication times between 
the nodes and the middle station over optical fibres, where Z¢ approx- 
imates the speed of light in an optical fibre. The factor > appears as, 
inthe implemented experimental sequences, the electronic delay for 
the atomic readout is only applied after a successful heralding event 
(Methods, Experimental sequence). This is not feasible when the two 
nodes are physically separated by a distance and, effectively, reduces 
the entanglement generation rate by a factor up to two for a physical 
separation L compared to the values observed in this work. 

For L <54km the entanglement generation rate is mainly reduced 
by the rapidly decreasing repetition rate R, for example, from L =O 
to L =33 km, R(L) drops by a factor of 8 while the success probability 
reduces by a factor 4.5. Only for distances L > 54 km will the exponential 
dependence of the success probability outweigh the dependence of 
the repetition rate. Extended Data Fig. 6 shows the expected entangle- 
ment generation rate according to equations (4)-(6) extrapolated to 
distances up to 100 km, together with the entanglement generation 
rate of the data presented in the main text. 

Note that the fidelities presented in the main text equal the expected 
fidelities for a physical separation of the nodes by a distanceL. Moreo- 
ver, the measurement times reported in the main text include a duty 
cycle of the experiment of roughly half, for example, the time effectively 
used for the entanglement generation tries approximates half of the 
reported measurement times (Methods, Experimental sequence). 


Atom-atom state fidelity 

The state of the atomic quantum memories is encoded in two magnetic 
sublevels of the rubidium ground state 5’S,/.|F = 1), which, however, is 
aspin-1 system. Besides the qubit states |m, = +1), also the state|m, = 0) 
can be populated, because of, for example, magnetic fields ina direc- 
tion not coinciding with the quantization axis. Hence, the atom-atom 
state effectively occupies a 3 x 3 state space. Assuming isotropic 
dephasing towards white noise, the fidelity relative to a maximally 
entangled state is therefore estimated as F> 1/9 + 8/9V , whereV is 
average visibility in three orthogonal bases. The commonly used fidel- 
ity estimation from visibilities for an entangled two-qubit system is 
F2=1/4+3/4V. For the effective two-qutrit system, however, this would 
result in an higher fidelity and would overestimate the fidelity of the 
generated atom-atom state. 


Polarization control of long fibres 

Stress or temperature induced polarization drifts in the long fibres are 
compensated using an automated polarization control. The polari- 
zation control is performed every 7 min, takes on average 20 s, and 
is based ona gradient descent optimization algorithm. In this way, 
polarization errors are kept below 1% during all measurements. 

The fibre polarization is optimized using laser light at the 
single-photon frequency with sufficient optical power to be detected 
by conventional photodiodes. Two polarization directions are used, 
vertical and diagonal linear polarizations, which are alternated at 10 Hz. 
The light is overlapped at the nodes with the complete single-photon 
path up to the detectors. In both output arms of the beamsplitter, a 
flip-mirror reflects the classical light into a polarimeter during the 
optimization. Three fibre polarization controllers are connected to 
the fibre beamsplitter of the BSM: at both input ports and at one out- 
put port, and are set according to the result of the gradient descent 
optimization algorithm. 

Polarization drifts in the long fibres are compensated in our setup, 
which includes a 700 metre fibre crossing public space anda four lane 
street. Recently, polarization drifts over a10 km field deployed fibre 
were characterized and compensated“. In our setup with a configura- 
tion of 32.4 km spooled and 0.7 km field fibre, we observe similar drifts. 
This indicates that a setup including longer field deployed fibres does 
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not introduce substantially more drifts than observed now and that 
the currently used system can compensate for polarization drifts in 
longer field deployed fibres. 


Entanglement swapping fidelity 

The single photons are detected with a BSM device consisting of a fibre 
beamsplitter, two polarizing beamsplitters (Wollaston Prisms) and 
four SNSPDs, as illustrated in Fig. 1 of the main text. In this setup, the 
fibre beamsplitter guarantees a unity spatial overlap of the photons 
originating from the nodes, whereas the polarizing beamsplitters and 
single-photon detectors allow for polarization analysis in both output 
ports. The detectors, labelled H,, V,, Hand V3, are not photon number 
resolving, and hence we only distinguish between six coincidence 
combinations of two detectors, see Extended Data Table 2. For the 
purpose of a BSM, we can categorize these combinations into three 
groups: D,, D_,and Dg. Here, detector combinations in group D, and 
D_herald the Bell states |*) and |W"), respectively, whereas combina- 
tions in group Dg should not occur for perfectly interfering photons 
and are discarded in the analysis. However, the relative occurrence of 
these events is used in the following to quantify the two-photon inter- 
ference contrast. 

For not interfering photons, two-photon events are evenly distrib- 
uted between the 16 possible detector combinations (not considering 
experimental imperfections). As the order of the detector combination 
is not of interest, for example, (V,,H,) is similar to (H,,V,), we end up with 
ten distinct coincidences and their probabilities, as listed in Extended 
Data Table 1. For perfectly interfering photons, the probabilities differ: 
the probability to detect the D, group vanishes and all four Bell states 
are detected with a probability of 1/4, whereby the |*) Bell states fall 
into the group ‘not detected’ for the setup used. 

The two-photon interference contrast is defined as” 


2Np, 


as ei 
Np, + Np_ 


(7) 


where WN, is the number of events in detection group k. With this defini- 
tion and the probabilities of the different coincidences, the contrast 
equals zero for not interfering photons and one for perfectly inter- 
fering photons. See ref. * for a thorough analysis of the two-photon 
interference contrast and entanglement swapping fidelity, including 
experimental imperfections. 

The interference contrast is measured as follows. During measure- 
mentrunsall single-photon detection events are recorded, which allows 
to count the number of occurrences of the coincidence events for all 
three detection groups. Next, the interference contrast is evaluated 
using equation (7). To verify this method, we additionally evaluate the 
contrast of not interfering photons. This is done by analysing coinci- 
dence detections of two photons originating from distinct entangle- 
ment generation tries. In this way, the photons did not interfere because 
the photon wave-packets are completely separated in time. Fig. 2b of 
the main text shows exactly this for the L = 6 km measurement. Shown 


are the normalized wrong coincidences, defined as 1-C, for varying 
time differences between the photons (Az). Note that the horizontal 
spacing of the measurement times equals the repetition rate of the 
entanglement generation tries. 

The entanglement swapping fidelity is mainly limited by two effects. 
First, it is limited by experimental imperfections that reduce the indis- 
tinguishability of the two photons, for example, as discussed in the main 
text, by an imperfect time overlap of the two photon wave-packets. 
Second, it is limited by double excitations stemming from the finite 
duration of the excitation pulse. For a detailed description, see refs.°**. 
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Extended Data Fig. 2 | Atom-photon entanglement distribution for Node1. 
The figues show the observed atom-photon state correlations for L, = 2.6, 5.4, 
11.3, and 16.5 km, as indicated in the labels at the top left corner of each 
subfigure. The measurements include an atomic readout delay to allow for 
two-way communication with the middle station. The error bars indicate 
statistical errors of one standard deviation. 
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Extended Data Fig. 3 | Atom-photon entanglement distribution for Node 2. 
The figues show the observed atom-photon state correlations for L, = 3.3, 11.4, 
and 16.6 km, as indicated in the labels at the top left corner of each subfigure. 
The measurements include an atomic readout delay to allow for two-way 
communication with the middle station. 
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Extended Data Fig. 4| Atomic state evolution in two bases. Atom-photon 
entanglement simulations (solid lines) and measurements (points) of Node1 
with photon detection at 780 nm using a5 mlong fibre. The atomic state 
readout orientation and time is varied to characterize the memory coherence. 
The labels in the legend indicate the|atom photon) state analysed. Simulation 
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parameters: trap waist 2.05ym, trap depth 2.32 mK, atom temperature 50K, 
bias field B, = 75.5 mG, and field fluctuations AB, = 0.5 mG Gaussian 
distributed. With these parameters, we observe de- and rephasing of the 
atomic state due to the longitudinal polarization components at the trap 
frequency of 70 kHzandalarmor precession at a frequency of 105 kHz. 
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Extended Data Fig. 5| Experimental sequence of the atom-atom 
entanglement generation. In both nodes asingle atom is trapped and cooled 
using polarization gradient cooling. Next, the entanglement generation tries 
start containing state preparation and synchronous excitation of the atoms. 
The atoms are re-cooled after 40 unsuccessful entanglement generation tries. 
The travel time of the photons to the middle station equals 14 to 84. us, 
depending onthe fibre length. With regard to the assumption that the 
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heralding signal is communicated back to the nodes along fibres of respective 
lengths, an additional delay is included to account for the communication time 
from the middle station back to the nodes. After a 200 ms interval of 
entanglement generation tries the presence of the atoms inthe trapsis 
checked using fluorescence collection with an APD (see text). The QFC takes 
place subsequent tothe atomic excitation. 
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Extended Data Fig. 6| Entanglement generation rate for varying link attenuation rate in optical fibres of 0.2 dB/km at 1517 nm (blue dashed line) and 
lengths. The plus markers show the observed entanglement generation rates 4.0 dB/km at 780 nm (gray dashed line). The difference between the blue solid 
(r=nR) for L = 6,11, 23, and 33 km, that is, the data presented in Fig. 3 and 4 of and dashed lines originates from the reduced repetition rate R(L) due tothe 
the main text. The blue solid line is asimulation of theentanglement generation longer communication times between the nodes and the middle station. For 
rate based on Equations 4-6 and extrapolated up toZ = 100 km. The dashed L<54km, the reduced repetition rate is the main contributor to the lower 


lines are proportional to the decay in success probability («n(L)) for an entanglement generation rate. 


Extended Data Table 1| Long fibre configurations and 
corresponding atomic readout times 


L(km) Li (km) La(km) Ai (JB) Az (GB) ti (us) ta (us) 
6 2.6 3.3 -0.7 0.8 28.5 35.5 
om 5.4 5.5 1.5 “1.3 57.1 71.0 
23 11.3 11.4 -3.3 -2.8 114.2 124.3 
33 16.5 16.6 45 44 171.2 177.5 


The fibre link lengths and corresponding atomic readout times for the experiments presented 
in the main text and the atom-photon characterization measurements. L, (L,) equals the 

fibre length between Node 1 (Node 2) and the middle station. A, (A,) gives the attenuation in 
the fibre network between the node and the middle station, this includes inefficiencies of 
fibre-to-fibre connectors. 
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Extended Data Table 2 | Possible two-photon coincidences 


detection coincidence nointerference __ perfect interference 
Hi,Hi 1/16 1/8 
H2,He 1/16 1/8 
not detected 
ViVi 1/16 1/8 
V2,V2 1/16 1/8 
Hi,H2 1/8 0 
Dg — discarded 
Vi,V2 1/8 0 
Hi.Vi 1/8 1/8 
Dy > |wt) 
H2,V2 1/8 1/8 
Hi,V2 1/8 1/8 
D_ >|) 
Vi,He 1/8 1/8 


The probabilities per coincidence event are given for not interfering photons and perfectly 
interfering photons. 


